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ABSTRACT 

We present theoretical calculations for the differential distribution of stellar orbital 
eccentricity for a sample of solar-neighbour halo stars. Two types of static, spherical 
gravitational potentials are adopted to define the eccentricity e for given energy E 
and angular momentum L, such as an isochrone potential and a Navarro-Frenk- White 
potential that can serve as two extreme ends covering in-between any realistic potential 
of the Milky Way halo. The solar-neighbour eccentricity distribution AN(e) is then 
formulated, based on a static distribution function of the form f(E,L) in which the 
velocity anisotropy parameter f3 monotonically increases in the radial direction away 
from the galaxy center, such that (3 is below unity (near isotropic velocity dispersion) 
in the central region and asymptotically approaches ~ 1 (radially anisotropic velocity 
dispersion) in the far distant region of the halo. We find that A7V(e) sensitively depends 
upon the radial profile of /3, and this sensitivity is used to constrain such profile 
in comparison with some observational properties of Ai\T b s (e) recently reported by 
ICarollo et al.l (|2010t ). Especially, the linear e-distribution and the fraction of higher-e 
stars for their sample of solar-neighbo ur inner- halo stars rule out a constant profile 
of f3, contrary to the opposite claim bv lBond et al.l ([2010). Our constraint of (3 < 0.5 
at the galaxy center indicates that the violent relaxation that has acted on the inner 
halo is effective within a scale radius of ~ 10 kpc from the galaxy center. We discuss 
that our result would help understand the formation and evolution of the Milky Way 
halo. 

Key words: methods: analytical - Galaxy: evolution - Galaxy: formation - Galaxy: 
halo - Galaxy: kinematics and dynamics. 



1 INTRODUCTION would provide us with a clue for looking into the early dy- 

_ , , . . „ T . . , . _ namical state of the Milky Way. 
Stellar halo of the Milky Way contains invaluable informa- 
tion to understand the Galactic formation and evolution, be- 
cause its collisionless nature has preserved some kinematic The first major milestone in the statistical study of or- 
properties of individual stars in the halo. The most com- bital eccentricity of halo stars was made by Eggen, Lynden- 
monly probed among them is the orbital eccentricity which Bell & Sandage (1962), who used a sample of 221 solar- 
remains almost unchanged as long as the gravitational po- neighbour stars with available (U, V, V^)-velocity data. They 
tential changes adiabatically. This quasi-adiabatic invari- calculated the model orbits of individual stars and found a 
ance of eccentricity leads to an idea that the eccentricity correlation between ultraviolet excess and orbital eccentric- 
distribution of halo stars has been conserved until present ity, such that metal-poor stars, which constitute the stellar 
since the end of the last rapid change of the potential. Such halo, are in eccentric orbits with complete lack of nearly 
rapid change may be related to some dynamical events - circular orbits. They interpreted this as the evidence that 
such as the major merger - in the course of forming the the stellar halo of the Milky Way formed out of a rapidly 
Milky Way and m ay have caused th e violent relaxation of collapsing proto-Galactic gas cloud. However, their sample 
the stellar system l|Lvnden-Bellll 1967t) . Thus, detailed anal- of halo stars, mostly taken from the catalog of high proper 
ysis of the observed eccentricity distribution of halo stars motions, turned out to be biased toward eccentric orbits. In 

fact, based on non-kinematically selected samples, various 
authors did show the existence of halo stars in nearly cir- 

* E-mail: khattori@ioa.s.u-tokyo. ac.jp (KH); yoshiiOioa.s.u- cular orbits an d required th e origi nal interpretation of the 

tokyo.ac.jp (YY) halo collapse bv lEggen et alj (|l962T ) to be significantly mod- 
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ified llYoshii fc Said Il979l ; Norris, Bessell, & Pickles 1985; 
IChiba fc Yoshiill 19971 ). 

Some attempts have been made along to explain the 
obse rved eccentricity distr ibution of halo stars. For exam- 
ple, IChiba fc Yoshiil |l998|) showed that a simple form of 
the distribution function that allows a radially-anisotropic 
velocity dispersion is more or less consistent with the ec- 
centricity distribution obtained from the Hipparcos data. 
However, because of the small sample size, they could not 
go further to fully interpret the eccentricity distribution, i.e., 
they could not well constrain the global dynamical structure 
over the halo by using the eccentric orbits that could travel 
far beyond the solar neighbourhood. 

Recently, with advent of huge surveys including the 
SDSS (Sloan Digital Sky Survey), much larger samples of 
halo stars are available, and it has become possible to discuss 
the observed eccen tricity distribution o n a firm statistical 
basis. In particular. ICarollo et al.l (|201dl ) analysed ~ 10, 000 
SDSS spectra of non-kinematically selected calibration stars, 
and derived their metallicities and radial velocities with rea- 
sonable accuracy. Using these as well as the proper motions 
from other catalogs, they found a marked metallicity depen- 
dence in the eccentricity distributions of halo stars obtained 
in different ranges of metallicity. They argued that this 
metallicity dependence would have arisen from a distance- 
dependent mixture ratio of two distinct stellar components 
of inner and outer halos having different eccentricity distri- 
butions and different metallicities. While their argument is 
intriguing in itself, more theoretical works are necessary to 
have insight into the origin of the dichotomic halo. 

As a first step toward impr oving current und erstanding 
of the eccentricity distribution. lHattori fc Yo shii (2010; HY 
hereafter) performed theoretical calculations of eccentricity 
distribution for the whole halo by modeling the gravitational 
potential and the distribution function. When applying to 
the solar-neighbour halo stars, however, the analysis in HY 
has to be modified, because such a sample is biased against 
stars that spend less fraction of orbital period in the solar 
neighbourhood. In this paper, by explicitly taking this bias 
into account, we re- formulate the analysis in HY and prop- 
erly compare our theoretical distribution to that observed 
for the solar-neighbour halo stars. In this way, we would 
be able to constrain the halo dynamics hopefully back to 
an epoch of violent relaxation that occurred in forming the 
Milky Way. 

In section 2, we present our new formulation for theo- 
retical eccentricity distribution of the solar- neighbour halo 
stars AN(e), with an emphasis on the difference from the 
formulation in HY. In section 3, we calculate AiV(e) for two 
spherical models of the Milky Way halo. In section 4, we 
demonstrate how our result in section 3 is beneficial to in- 
terpret the solar-neighbour data. In section 5, we summarize 
the basic ideas from our analysis. 



2 FORMULATION 
2.1 Survey region 

As in HY, we assume that the stellar halo is a spherical 
system that consists of many halo stars in a steady, spherical 
gravitational potential of the dark matter halo. We further 



assume that kinematic data of halo stars can be obtained 
only in a shell-like region defined as 

r_ < r < r+, (1) 

where r = \r\ is the Galactocentric distance. We shall here- 
after call this the survey region. We note that the actual sur- 
vey region in this shell should be confined in the solar neigh- 
bourhood, centered at the position of the sun. However, this 
realistic definition of the actual survey region barely makes 
any difference to any result from the shell-like region [in- 
equality (JTJ] as far as the halo is assumed to be of spherical 
symmetry in our analysis. See also AppendixfX]and footnote 

EE! 

2.2 Stellar orbital eccentricity in a model halo 

Given a spherical halo potential V(r), the orbital eccentric- 
ity of a star is practically defined as 

r^TDO ^*Dcri 



(2) 

^apo "T ^pcri 

where r apo and r por i are the apocenter and pericenter dis- 
tances, respectively, and are given by two real solutions 
(7"a P o > r pal i) of the following equation: 



E = V(r) + 



id 

2r 2 



V eS (L;r), 



(3) 



where E and L are the specific energy and angular momen- 
tum, respectively. As discussed in Appendix A of HY, bound 
stars (E < 0), provided they are observable in the survey 
region, reside only in a limited region of the (E, L)-phase 
space, and we shall concentrate on these stars below. 



2.3 Differential eccentricity distribution of 
observable stars 

The distribution function for a system of spherical symme- 
try takes the form f(E, L) ] according to the s trong Jeans 
theorem (|Lvnden-Belll 1 19601 ; iLvnden-Belll I1962T I . By chang- 
ing variables and integrating over spherical coordinates, we 
obtain the number of stars in the survey region [inequal- 
ity fTJ] and m a phase space volume dEdL 2 at an epoch of 
observation: 

AN(E, L) dEdL 2 = An 2 f{E, L)AT r (E, L) dEdL 2 . (4) 

This equation is essentially the same as equation (7) of HY, 
except that the radial period T r in HY is replaced by the 
observable time AT r which is the total length of time per 
radial period a star with a given E and L spends inside the 
survey region: 

AT r (E, L) = 

/■miii(r + ,r apo ) 

2#(min(r+, r apo ) - max(r_, r pcri )) / -, — r, (5) 

imax(r_ ,r pcri ) \ V r\ 

where 9{x) is Heaviside step function. Note that observable 
stars that satisfy 



r- < r-apo and j- pcri < rj, 



(6) 



pass through the survey region at an epoch of observation 
with a probability AT r /T r . In particular, observable stars 
with 



r~- <C ?"p er i *-C /"apo ^ r*-(- 



(7) 
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are always inside the survey region and AT r coincides with 
T r . On the other hand, unobservable stars, for which 

r apo < r- or r+ < r pcri , (8) 

do not pass through the survey region and therefore AT r = 
0, as clearly seen from equation (J5]). 

By using these quantities and remembering that L 2 is a 
function of E and e, we obtain the independent differential 
eccentricity distribution of halo stars in the survey region 
[inequality fT}]: 



An(E,e) = An 2 f(E, L)AT r (E, L) 



de 



(9) 



This expression differs from HY, i.e., we have included the 
distribution function f(E,L), instead of placing it outside 
as a weight to be added to the _E-dependent eccentricity 
distribution in equation (8) of HY. We then express the dif- 
ferential eccentricity distribution as 



AJV(e) = 



An(E,e) dE. 



(10) 



V(r_) 



Here, the integral over E should be performed over the en- 
ergy range of bound, observable stars with AT r > 0. The 
allowed iJ-region for such stars is presented in the (E, L) 
diagram in Appendix 1X1 



In this paper, following ICuddefordl | 1991), we model the 
stellar halo as a family of distribution functions of the formQ 



f(E,L) = 



\g{Q)L 
l°> 



-200 



if (E, L) is 'allowed' 
otherwise, 



(12) 



where g(Q) is an arbitrary function of Q which is a linear 
combination of E and L 2 [see section [3.1.31 below] . and Po is 
P at the galaxy center. In this model, the radial profile of P 
as a function of r is expressed as 



P = Pa 



+ 



t + r„ 



(13) 



where r a is the so-called anisotropy radius which character- 
izes the profile of p. [In section[3l we re-arrange this form in 
terms of Po and /?© instead of Po and r a -] 

The profile of p for r a = oo reduces to a constant profile 
with /3 — Po- On the other hand, the profile for a finite 
value of r a behaves such that f3 increases monotonically from 
ft at r = to /J = 1 at r > r a . The functional form 
of the profile in equation (|13[l . though not covering /3 that 
decreases with increasing r, is enough to examine effects of 
the radial gradient in P upon An(E,e) as well as AJV(e). 
In a forthcoming paper (Hattori 2011, in preparation), we 
will discuss another family of distribution functions in which 
Po > Pq is also allowed. 



2.4 Distribution function of stellar halo 

Our knowledge of the distribution function of the stellar 
halo f(E, L) is fragmentary. Observationally, it is indicated 
that the first moments or the mean values of velocity com- 
pone nts are near zero ex cept for the rotation component 
(e.g. IChiba fe"l3 ccrs 2000). [In fact, the first moments for 
the spherical model with no systematic rotation are zero.] 
Thus, the most useful constraint on f(E, L) is the behaviour 
of the velocity anisotropy which is usually parameterised as 



1-^1 

9 ' 



(11) 



where ay is the radial velocity dispersion and a t is the tan- 
gential velocity dispersion projected onto the spherical 8-(j> 
surface. 

Estimates of p in the solar neighbourhood are rather 
conve rgent and are in a range of /3q = 0.4 — 0.7 
(e.g.. lYoshii fc Saidll979l ; Gilmore, Wyse & Kuijken 1989; 
ICarollo et alj|20ich. while its value beyond the solar circle 
is less certain. ISommer-Larsen et al.l (|l997f l used a sample 
of field horizontal branch stars and claimed that /3 decreases 
with increasing r (i.e., the tangential anisotropy in velocity 
dispersions dominates in the far distant halo), which they 
regarded as a sign for the bot tom-up galaxy formation sce- 
nario, while iBond et all ||2010T ) used a stellar sample of SDSS 
and suggested that /3 is more or less constant within 5 kpc 
of the sun. On the other hand, iV-body simulations suggest 
that the cold collapse of the halo, which is most likely to 
trigger the violen t relaxation, favours P that inc reases with 
increasing r (e.g. Ivan Alb ada 1982; Voglis 1994). Similarly, 
recent iV-body + gas-dynamical simulations suggest that P 
of a stellar halo increases with increasing r if it is formed 
through disruption of satellite galaxies, (e.g. Abadi, Navarro 
& Steinmetz 2006). 



2.5 Limiting cases of the survey region 

As mentioned above, the finite radial thickness jr + — r_ of 
the survey region affects directly AN(e), but the difference 
shows up only through AT r . To see further this effect, it is 
instructive to consider two limiting cases: (1) the large-shell 
limit of r_ = and r+ — > oo, and (2) the thin-shell limit of 
r_ = r-0 — 5r/2 and r+ = r© + 5r/2 with 8r — > 0, where r© 
is the radial distance to an observer from the galaxy center. 

Obviously, in the case (1), it follows that AT r — T r 
and the result of AN(e) is the same as N(e) given in HY. 
Therefore, in this subsection, we do not repeat the result of 
this case. In the case (2) of thin-shell limit, use of equation 
(O gives 

St 

AT r -> 2 0(r apo - r )(9(r Q - r por i), (14) 

VrQ 

where v r Q = \v r \ is the radial speed at r = r© which de- 
pends on E and e0 In the following discussion, we consider 
the eccentricity distribution, using a distribution function 
of the above form [equation l |12p ] with a constant profile of 
P (r a = oo) in two extreme gravitational potentials of the 
central point mass and the truncated homogeneous density 
distribution. We note that these potentials allow analytic 
expression of Anp(E,e) and serve as useful reference. Here 
and hereafter, the subscript /3 to An(E, e) or AJV(e) stands 
for the case of constant profile of p. For our explanation be- 
low, it is useful to introduce the specific energy of a star in 



1 The allowed region of the (E, L)-phase space for bound stars 
is described in Appendix A of HY. 

2 Strictly speaking, equation H14I I does not hold when v t q = 0. 
However, we neglect such a physically trivial exception in our 
analysis. 
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a circular orbit with orbital radius r: 
E c (r) = V(r) + ±r±V(r). 



2. 5. 1 Central point mass 



(15) 



The gravitational potential arising from the central point 
mass is Keplerian and is given by V(r) = —GM/r, where 
M is the mass and G is the gravitational constant. By eval- 
uating r apo and r pcr i, we obtain 



r o)9(ro - r P eri) = 



E 
Eq~ 



1 



(16) 



where Eq = E c (r Q ) = -GM/(2r Q ). 

When E < V(r@), equation (|16|l vanishes, therefore 
AT r = and Anp — 0. On the other hand, when V(tq) < 
E < 0, AT r can be expressed as 



AT r 



ASr 0(e - e cu t) 



(17) 



which allows the analytic expression of Anp{E, e) for a con- 
stant profile of ft as follows: 



An fj (E,e) 



32w 2 5r(GM) 



2-2/3 



9(E) 



?(e — e cu t) 



(18) 



where the cutoff eccentricity e cu t is defined as 

e cut = \E/E @ - 1|. (19) 

Evidently, e cut gives the same value for a pair of two energies 
E and W which are related to each other via E — Eq = Eq — 
W. Therefore, we see from equation (|18[) that the shape of 
Anp(E, e) for a given E is identical with that of Anp(W, e) 
for the corresponding W . Because of this relation, we only 
need to know the shape of Anp(E, e) for either E Eq or 
E ^ Eq. Results of An@(E, e) for several combinations of /3 
and E are shown on the left-hand panels of Figure [T] 

2.5.2 Truncated homogeneous sphere 

A homogeneous density distribution within truncated sphere 
is expressed as p(r) = 3M/(47rr 3 ) at r < r t and p(r) = 
at r > r t , where M is the total mass of dark halo and r t 
is the truncation radius. The gravitational potential arising 
from this density distribution is given by V(r) — E m i n + 
GMr 2 /(2rl) at r < r t , where E min = -3GM/(2r t ). If the 
stellar system is confined inside the truncation radius r t , by 
evaluating r apo and r pcr i, we obtain 



^(rapo - r Q )8(rQ - r pcli ) = i 



2e 



1 + e 2 



\E-E Cz 

E — Em[ 



(20) 



where E e = E c (r e ) = E min + GMr%/rt 

When E < V(ro), equation (|2U|) vanishes, and therefore 
AT r = and An^ = 0. On the other hand, when V(ro) < 
E < V(rt), equation (|20p reduces to 8(e — e cut ) and thus 
AT r can be expressed as 



AT r 



25r 



y/GM/n Vq 



E — Eg 



E — Err 



2e 



- , 



(21) 



which allows the analytic expression of Anp(E, e) for a con- 
stant profile of j3 as follows: 



Anp(E,e) = 

2 5rf rl \ 3/2 " /3 

e(l-e 



{E @ - E min ){E - E min f-^g{E) 

2x1-2/3 



0(e - ecut) 



where the cutoff eccentricity e cut is defined as 



f, (22) 



e cu t 



E — Emi 

\E-Er. 



E — E m j 
E-Ec. 



(23) 



Similarly to the case of central point mass, e cut gives the 
same value for a pair of two energies E and W which are 
related to each other via (E — Emi n )/(E — Eq) — (W — 
-Emin) / (Eq — W). Therefore, we see from equation (|22p that 
the shape of Anp(E, e) for a given E is identical with that of 
Anp (W, e) for the corresponding W . Because of this relation, 
we only need to know the shape of Anp(E,e) for either 
E ^ Eq or E ^ Eq, although we have to keep in mind 
that we truncate the stellar system at E = V(r t ). Results 
of Anp(E, e) with r© = 0.15r t for several combinations of j3 
and E are shown on the right-hand panels of Figure [T] 

2.5.3 Insights from the thin-shell limit 

In HY, ng(E, e) for the whole stellar system with a constant 
profile of /3 is found to be given by F(E) x G(e) in the 
extreme gravitational potentials considered here, so that the 
shape of np(E,e) does not depend on E (see Figure 1 of 
HY). On the other hand, Anp(E,e) for the survey region 
described above depends on E. Since Anp/np tx AT r /T r [f] 
and T r does not depend on eQ extra e-dependence of Anp 
over np is attributed exclusively to the e-dependence of AT r . 

In the thin-shell limit, AT r shows a marked e- 
dependence, especially near e = e cut : it vanishes at ^ 
e < e C ut, diverges at e = e cut + 0, and rapidly decreases at 
e > e cu t, with its slope becoming less and less steep toward 
higher e. This marked e-dependence of AT r is reflected to 
Anp and the shape of Anp at e < e cut is essentially deter- 
mined by AT r , regardless of /?. At higher-e region, however, 
the L _2/3 term in the distribution function [equation (|12[) ] 
starts to affect Anp with its slope continuing to decrease 
toward higher e for j3 < 0.5, while showing upturn there for 
P £ 0.5. 

In section 2.4 of HY, we have shown that np for more 
centrally concentrated system is more weighted toward high- 
e region. This trend can be confirmed for Anp shown in Fig- 
ure[T] We see from this figure that for given /3 and E, Anp for 
the central point-mass model is more weighted toward high- 
e region than that for the truncated homogeneous sphere 

3 Notice that a way of defining np without g(E) as in equation 
(8) of HY differs from Anp which includes g(Q) as in equations 
||9} of this paper. 

4 This is the case for these extreme potentials considered here as 
well as in the isochrone potential below [see equations (13), (23), 
and (33) of HY], In general, T r depends on both E and e, but 
e-dependence is often very weak. 
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model. However, Figure [T] also suggests that this trend is 
not very prominent for a stellar system with ft ~ 1. 

From the analytic expression of e cu t, we see that Anp 
depends not only on specific energy E of observed star^f] but 
also on Eq which is a function of M and vq for the case of 
central point mass, or a function of M, r Q , and r t for the 
case of truncated homogeneous sphere. Thus, realistic mod- 
eling of the Milky Way is necessary to perform theoretical 
calculations of An^ expected at the position of the sun. 



3 ECCENTRICITY DISTRIBUTION OF 
SOLAR-NEIGHBOUR HALO STARS 

The formulation in the previous section can be applied to 
more general gravitational potentials, such as the isochrone 
and NFW potentials, which we consider to be limiting cases 
of reasonable descriptions of the mass distribution in dark 
haloes. In this section, we derive AiV(e) for these potentials 
by specifying the survey region and the distribution func- 
tion. 



3.1 Calculation setup 

3.1.1 Gravitational potential 

The gravitational potential of the isochrone model (|Henonl 
1 19591 ) is given by V(r) = —GM/(b + Vb 2 + r 2 ), where M 
is the total mass and b is the scale length parameter . For 
the Milky Way halo, following Ichiba fc YoshiH l|l998j ). we 
set b = 5.2 kpc and s /GM/{2b) = 385 km s" 1 . 

The NFW density profile of dark matter produces a 
gravitational potential of the form V(r) = — AirGpoa 3 [In (1 + 
r/a)]/r, where po is the typical density and a is the scale 
length parameter (Navarro, Frenk & White 1997). Consider- 
ing typical values for the Milky Way-size dark matter haloe s 
in recent cosmological simulations (e.g., IZhao et al.l 12003). 
we set a = 20 kpc and <*J A-KGpoa 2 = 360 km s _1 , together 
with additional values of r2oo = 200 kpc for the virial radius 
and c = r2oo/a = 10 for the concentration parameter. 



3.1.2 Survey region 

The survey region adopted is given by r_ = 7 kpc and 
r+ = 10 kpc in inequality {TJ. This region matches with 
the sp atial criteria for sampling halo stars bv lCarollo et al.l 
l|2010| ). namely, 7 kpc < R < 10 kpc (R is the projected 
Galactocentric distance onto the disk plane) and d < 4 kpc 
(d is the distance from the sun). 



3.1.3 Distribution function 

The distribution function adopted is given by the form 
in equation (|12[1 . with radial profile of /3 which is re- 
parameterised in terms of /3o at the center and /3q at rg, 
instead of f3o and r a in equation (|13[1 . In our analysis be- 
low, we consider several combinations of (/3o,/3©) provided 
^ /So ^ /3q. In particular, for fio = /9q, a constant profile 
results. 



For the independence of e cu t, see Appendix [A] 



As for g(Q) with the argument defined as 
L 2 /(2r 2 ), we adopt 



g(Q) 



A (exp | 
0, 



2-S.l 



-1 



if Q < E t 
otherwise, 



E + 



(24) 



where A is a constant, Et is a truncation energy, a is a 
velocity parameter described below. 

In general, truncation of the distribution function at 
Q — E t for a stellar system in a spherical potential V(r) 
guarantees stars to be confined within a truncation ra- 
dius rt for which V(rt) = -EtQ For the isochrone po- 
tential, we express the dimensionless truncation energy as 
e t = 2bEt/(GM). Using r t = 200 kpc as a characteris- 
tic value of reference, we obtain e t = —0.05. On the other 
hand, for the NFW potential, as in HY, we express it as 
e t = E t /{4:nGpQa 2 ) = — [ln(l + c)]/c, by identifying the 
virial radius r2oo with the boundary of dark matter halo. 
Then, using c = 10 (section l3.1.ip . we obtain e t = —0.23. 

In our model, the radial velocity dispersion ay decreases 
with increasing r, which is consistent with some observations 
(e.g. iBrown et al.l |2010| ) . In addition, under a fixed set of 
the gravitational potential and the truncation energy Et, 
the radial profile of a r /a is determined by /3q onljjj and 
this /?o-dependence is fairly weak(f| In this paper, a is set to 
be 200 km s _1 so that oy decreases with r from the central 
region, reaches ~ 160 km s" 1 in the solar neighborhood, and 
vanishes at r = r t . 



3.2 Results 

3.2.1 Energy-dependent eccentricity distribution An(E,e) 
for the stellar halo 

The central part of deriving An(E, e) lies in the observable 
time AT r . For the isochrone potential, after tedious algebra, 
we obtain the analytic expression of AT r as follows: 



AT r 



GM 



ex 2 + 2V1 + X 2 + 



min(r apo ,r + ) max(r pcri , r_) 



(l-e 2 ) 2 -(l + e 2 ) v / (l- e ^) 2 +4^ e 2 



1/2 



+ - 



1 ( Vl+^-V^-y'H^ 



-. x— min(r apo ,r + )/6 



(-£) 



3/2 



J x— max(r pcr i ,*r_ )/b 



(25) 



e This is because r apo of a star increases with E when L 2 is 
fixed, while it decreases with L 2 when E is fixed. 

7 The radial profile of at I a is determined by r a as well as 0o, 
since a 2 /a 2 = (1 - 0o)r 2 J(r 2 a + r 2 ). 

8 For -1/2 < A) < 1/2, cr r /(j is different from that for /3 = by 
only 8% (relative difference). This percentage becomes no larger 
than 19% for -1< fa < 1. 
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- Point-mass model 

- P = -1 
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Figure 1. Energy-dependent differential distribution An^(B,e) of stellar orbital eccentricity in a spherical system with constant radial 
/3-profile for an observer at tq = 8.5 kpc whose survey region is limited to nearby stars only (thin-shell limit). Shown are the results 
in two extreme models of gravitational potential, such as the point-mass model for several values of /3 on the left-hand panels (section 
I2.5.1|l and the homogeneous model with r^/rt = 0.15 for several values of /3 on the right-hand panels ( section I2.5.2J I, Shown by lines 
on each panel for given /3 are the results for several values of specific energy E equated to E c (r) of a star in circular orbit at different 
orbital radius r. Note that An^(B,e) is normalized such that Anp(E,e)de = 1. 



where e = 2bE/(GM), x pcli = r pcli /b, and x apo = r apo /b|f| 
Analytic expressions of L 2 and (dL 2 /de) E in terms of s and 

9 Note that when r_ < r pcr i < r apo < r+ is satisfied, the 
first term in the square bracket of equation l|25| l vanishes and the 
second term is equal to 7r/(— e) 3 ^ 2 , so that it gives AT r = T r = 
2nGM/{-2E) 3/2 [cf. equation (33) of HY]. 



e, which are necessary to derive An(E, e), are found in equa- 
tions (34) and (35) of HY. By substituting these quantities 
into equation (|9]) and using the parameters given in section 
13.1.11 we obtain the fully analytic expression of An(E,e). 



Alternatively, for the NFW potential, we obtain the for- 
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mal expression of AT r as follows: 



AT r = 



r = \J 27r G po e ( mill ( r +- r »p») ~ max(r_, r peri )) 



L 



min(r + ,r apo )/a 



xdx 



*(*•_, r peri )/a £x 2 +xln(l + l) 



(26) 



where e = _E/(47rGpoa 2 ) and A = L 2 /(4nGpoa 4 ). Analytic 
expressions of other necessary quantities are found in equa- 
tions (49) and (50) of HY, and we numerically calculate 
An(E,e). 

Figure [5] shows calculations of An@(E,e) for constant 
profiles of f3 = 0.4, 0.5, and 0.7. Three panels of differ- 
ent values of /? in the left-hand column show the results 
for the isochrone potential, and three panels in the right- 
hand column for the NFW potential. In different panels 
for different values of /3, shown by lines are the results 
for E = 0.98V(r_), E c (r Q ), E c (2r Q ), and E c (4r ), where 
E c (r) is defined in equation (|15|) . Irrespective of the poten- 
tial adopted, we see common features in Anp(E,e), i.e., its 
shape is weakly dependent on P in the low-e region, while the 
fraction of stars with high eccentricity notably increases as 
P increases. We see from Figure [2] that, just as in the case of 
thin-shell limit (Figure [TJ, Anp has a non-zero cutoff eccen- 
tricity e cu t below which Ana vanishes, when E < E c (r~) 
and E > E c (r+). This cutoff corresponds to r apo = r_ 
when E < _E c (r_), and r pcr i = r+ when E > E c (r+). On 
the other hand, we also see that there are some eccentric- 
ities at which Anp shows non-smooth behaviour. This sin- 
gularity of Anp corresponds to r apo = r+ or r por i = r_ 
and is analogous to the divergence of Ana at e = e cut 
in the thin-shell limit that corresponds to r© = r por i or 
r^po- The singularity associated with r apo = r+ appears 
when V(r + ) < E < E c (r + ) [e.g., E — _E c (r Q ) in our re- 
sults], while that associated with r pcr i = r_ appears when 
£ c (r_) < E [e.g., E = E c (r e ), E c (2r Q ), and E c (4r e )]. Both 
of these singularities appear in the overlapping _E-region of 
m^{E c (r-),V(r+)} < E < E c (r+) [e.g., E = £ c (r )B 
while none of them appears when E < min{_E c (r_), V(r+)} 
[e.g., E = 0.98V(r_)]. 



Particularly, as long as /3 — 0.4 — 0.7, AiV(e) is nearly linear 
in e in the low-e region (e < 0.5). 

Second, Figure [4] shows AJV(e) with non-constant pro- 
file of P for several combinations of (Pq, Pq). The linear trend 
of AJV(e) in the low-e region, which is seen in Figure |3j re- 
mains as long as /3© = 0.4 — 0.7 and Pq > 0. Figure [4] also 
shows that the fraction of high-e stars increases as Pq in- 
creases from to Pq. In particular, the fraction of stars 
with 0.9 < e < 1 for pa = Pq is nearly 2-3 times larger than 
that for p = 0. 

Third, Figure [5] shows AJV(e) made up of two 
stellar components. Following the recent observation by 
ICarollo et all (|2010h that the Milky Way halo consists of 
at least two components, such as the inner halo with 
Pq ~ 0.7 and the outer halo with Pq ~ 0.4, we here con- 
sider three cases: (i) a mixture of an 80% component with 
(/3 ,/3 Q )=(0.7,0.7) and a 20% component with (0.0,0.4), (ii) 
a mixture of an 80% component with (0.7,0.7) and a 20% 
component with (0.4,0.4), and (iii) a single component with 
(0.7,0.7). Calculating the eccentricity distributions for indi- 
vidual components, and adding them up according to the 
mixture ratio adopted, we obtain AN(e) for each of the 
cases (i) and (ii), to be compared with the case (iii). We see 
from this figure that contamination of the lower-/?© outer- 
halo component by 20% level gives no significant effect on 
AiV(e) for the higher-/?© inner-halo component, particularly, 
at e > 0.7. In other words, from the observed AiV(e) we can 
best constrain the profile of P for the inner-halo component, 
in a manner largely unaffected by possible contamination of 
the outer-halo component. 



4 DISCUSSION 

In the previous section, we have performed theoretical cal- 
culations of AJV(e) for the stellar halo, assuming spherical 
symmetry for simplicity. In this section, we demonstrate how 
our calculations would be useful in interpreting the observa- 
tions of kinematics of halo stars in the Milky Way. 



3.2.2 Eccentricity distribution AN(e) for the stellar halo 

Integrating An(E, e) over E [equation (fT0)l ]. we obtain 
AiV(e) for the isochrone potential and the NFW potential. 
The results of AN(e) for both potentials similarly show a 
notable dependence on the profile of p. In the following, we 
summarize the characteristics of AiV(e), which are common 
in both potentials except for some minor differences. 

First, Figure [3] shows AiVs(e) with constant profile of P 
for several values of p. A strong dependence of ANp(e) on P 
is apparent; As p increases, A/Vg(e) changes from a hump- 
like distribution (P < 0.4) through a nearly linear distribu- 
tion (P ~ 0.5 — 0.6), then to a steeper turnup distribution 
with a prominent increase in the high-e region (/? > 0.7). 



10 Anp(E c (rQ), e) has two singular points at e ~ 0.2 for both 
potential models. However, for the case of isochrone model, occa- 
sional proximity of these two singular points makes it difficult to 
see them separately in Figure [2] 



4.1 Comparison with the observed eccentricity 
distribution 

4-1.1 Summary of the latest observation 

So far the largest, kinem atically unbia s ed sa mple of halo 
stars is that obtained by ICarollo et al.l (|20 10T ) . They eval- 
uated eccentricities of solar-neighbour halo stars sampled 
in 7 kpc < R < 10 kpc and d < 4 kpc, using an ax- 
isymmetric Stackel-type poten tial which was first used by 
ISommer-Larsen &: Zhenl (|l990h for the analysis of stellar 
halo. Their finding is that the eccentricity distribution of 
stars differs significantly depending on the metallicity range 
chosen (see their Figure 5). The results for likely-halo stars 
with [Fe/H]< —1.5 and \z\ > 1 kpc (\z\ is the distance from 
the disk plane) are summarized as follows: (1) AA r b s (e) for 
stars with —2.0 <[Fe/H]< —1.5 shows a linear trend over 
the full range of eccentricity from e = to 1. (2) AA r b s (e) 
for stars with [Fe/H]< —2.0 also shows a linear trend at 
e < 0.6 — 0.8, but with a shallower slope toward higher e, in 
contrast with the case of -2.0 <[Fe/H]< -1.5. (3) AAT obs (e) 
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Figure 2. Energy-dependent differential distribution An^(B,e) of stellar orbital eccentricity in a spherical system with constant radial 
/3-profile for an observer whose survey region is 7 kpc < r < 10 kpc. Shown are the results in two realistic models of the gravitational 
potential for the Milky Way halo, such as the isochrone model on the left-hand panels and the NFW model in the right-hand panels. 
Shown by lines on each panel for given f) are the results for several values of specific energy E equated to E c (r) of a star in circular orbit 
at different orbital radius r. Note that An^E, e) is normalized such that An^(_B, e)de = 1. 



as above is unaffected by the spatial criterion for sampling 
stars at either 1 kpc < \z\ < 2 kpc or 2 kpc < \z\ < 4 kpc. 

ICarollo et alj (|201Cf ) expected that this metallicity- 
dependent AjV b a (e) could be explained by the metallicity 
dependence of mixture ratio of the inner- and the outer-halo 
components. According to their decomposition analysis, the 
sample for -2.0 <[Fc/H]< -1.5 is mainly (> 80%) con- 
tributed by the inner-halo component. Although they did 



not perform a similar decomposition analysis on the sam- 
ple with [Fe/H]< —2.0, they supposed a smaller fraction 
of high-e stars in this very metal-poor sample than that 
for —2.0 <[Fc/H]< —1.5, because of larger contamination of 
outer-halo component which is expected to possess a smaller 
fraction of high-e stars than the inner-halo component. 

By comparing all these latest results of AiV bs(e) with 
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Figure 3. Differential distribution ANp(e) of stellar orbital eccentricity in a spherical system with constant radial /3-profile for an 
observer whose survey region is 7 kpc < r < 10 kpc, embedded in the isochrone potential (left-hand panel) and the NFW potential 
(right-hand panel). Shown by lines on each panel are the results for /3 = 0, 0.4, 0.55, 0.7, and 0.9, in order from top line to bottom line 
in the smaller-e region. Note that ANp(e) is normalized such that J ANp(e)de = 1. 



our theoretical calculations of AN(e), we discuss their im- 
plications in the global kinematical structure below. 

4-1.2 Linear trend in the low- eccentricity region 

A linear e-distribution obtained bv lCarollo et al.l (|2010h pre- 
vails at least at e < 0.6, which is common to both samples of 
-2 <[Fc/H]< -1.5 and [Fe/H]< -2.0 [results (1) and (2) in 
section [4. 1 . 1 J . Combined with their estimation offlm ~ 0.7 
for the inner halo and Pq ~ 0.4 for the outer halo^jthe ob- 
served linear trend in the low-e region is well reproduced 
by our theoretical e-distribution AiV(e) for both of the 
isochrone and NFW potentials, provided Po > 0. We expect 
that this result holds for any realistic spherical mass distri- 
bution models of the Milky Way halo, because the isochrone 
and NFW potentials form two extreme ends of such models. 
In addition, it may as well hold for non-spherical mass distri- 
bution models, because the equi-potential surface is rounder 
than the equi-density surface and the estimation of eccen- 
tricity is hardly affected by the non-sphericity of the po- 
tential at least in the low-e region. Thus, regardless of the 
mixture ratio of the inner halo (Pq ~ 0.7) and outer halo 
(/3q ~ 0.4) and regardless of either the spherical or non- 
spherical mass distribution assumed for the halo, we still 
expect that a nearly linear e-distribution would occur at 
e< 0.5 -0.6. 

4-1.3 A constraint on the radial /3 -profile for the inner 
halo 

Figure 5 of ICarollo et~ai1 (|2010h shows that the linear trend 
of AJV bs(e) extends beyond the low-e region up to e = 1 
for a sample of stars with —2.0 <[Fe/H]< —1.5 [result (1) in 
section [4.1.1| . This trend can be used to constrain the radial 

11 These values of /3q are derived from (ay , ag , <t$ ) = 
(160,102,83) kms" 1 for the inner halo, and (178, 149, 127) 
km s - 1 for the outer halo, which were reported in ICarollo et al.l 

toid ). 



profile of /3 for the inner halo which dominantly contributes 
to such a sample. 

In ge neral, when a non-s pherical potential is adopted, 
just as in ICarollo et ail (|201fjh . the estimation of eccentric- 
ity of a star tends to be systematically larger than that for 
a spherical potential. This discrepancy cannot be ignored 
at high-e region, so that adopting a non-spherical poten- 
tial instead of a spherical one would enhance the fraction of 
high-e stars, while suppressing the fraction of low-e stars. 
Accordingly, the approximately linear shape of AA^ b s (e) 
up to e = 1 for a sample of likely inner-halo sta r s wit h 
-2.0 <[Fe/H]< -1.5, reported by ICarollo et all |2010l ). 
would be changed into a hump-like shape when the eccen- 
tricity in their analysis is instead defined in our spherical 
model of the Milky Way halo. We then constrain the radial 
profile of P for the inner halo by examining how our spherical 
model could reproduce such a hump- like shape of AiV b s (e). 

Noting Pq ~ 0.7 for the inner halo, our models of 
AiV(e) for various /3-profiles with Pq = 0.7, shown on 
the bottom-row panels of Figure 0] would help constrain 
the /3-profile for the inner halo. We see from these pan- 
els that our models for more or less constant /3-profiles 
with Po — 0.5 — 0.7 and Pq = 0.7 would be ruled out. 
On the other hand, our models for notably r-dependent 
/3-profiles with Po < 0.5 and Pq = 0.7 would repro- 
duce AiV bs(e) if a non-spherical potential is assumed to 
be consistent with Carollo et al.'s analysis. Remembering 
(r a /r ) 2 = (1 — Pq)/(Pq — Po), a constraint of Po < 0.5 for 
the inner halo is equivalent to r a < 10 kpc. 

The key idea used a bove is that the samp le of stars with 
-2.0 <[Fe/H]< -1.5 bv lCarollo et all (|20ld ) has a smaller 
fraction of high-e stars when compared with that expected 
from our models with a constant profile of Pq = 0.7. Since 
smaller Pq results in a smaller fraction of high-e stars (see 
Figure [J), some contamination of a stellar component with 
smaller Pq would obviously decrease the fraction of high-e 
stars in the sample. Thus, we examine whether the observed 
small fraction of high-e stars can be explained by our mod- 
els if the contamination of the outer halo with Pq ~ 0.4 is 
taken into account. Figure [5] clearly shows this effect. For 
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Figure 4. Differential distribution AA^(e) of stellar orbital eccentricity in a spherical system with the /3- profile expressed by equation 
(1 1 3 1 > for an observer whose survey region is 7 kpc < r < 10 kpc, embedded in the isochrone potential (left-hand panel) and the NFW 
potential (right-hand panel). Shown by dot-connecting lines on each panel for given /3q are the results for several values of /3o- The case 
of /3q = /3q corresponds to the constant profile of 0. In particular, the case of 0q = 0.5 and /3© = 0.7 gives the nearly linear e-distribution 
over a full range of e for both potentials adopted. Note that AN(e) is normalized such that AN(e)de = 1. 



the cases (i) and (ii) in section [3.2.21 where 80% of the sam- 
ple is contributed by the inner halo with the constant pro- 
file of /3q — 0.7 and the other 20% by the outer halo with 
/3q = 0.4, the fraction of high-e stars in the sample becomes 
distinguishably smaller than that for the case (iii) for which 
the sample genuinely consists of the inner halo. However, 
we see from this figure that 20% contamination of the outer 
halo is not able to make a hump-like shape of AJV(e) and is 



not enough to explain the observed small fraction of high-e 
stars in the sample. Ther e fore, the observed e-distribution 
reported bv lCarollo et al.l (|2010h implies that the radial /3- 
profile for the inner halo is not constant. Rather, it is con- 
sistent with a /3-profile that increases away from the galaxy 
center with /3q < 0.5. 
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Figure 5. Differential distribution ANa(e) of stellar orbital eccentricity for stars in 7 kpc < r < 10 kpc that are made up of two stellar 
components, embedded in the isochrone potential (left-hand panel) and the NFW potential (right-hand panel). Shown are the results for 
three cases such as (i) a mixture of an 80% component with (/3n, /3©)=(0.7,0.7) and a 20% component with (0.0,0.4), (ii) a mixture of 
an 80% component with (0.7,0.7) and a 20% component with (0.4,0.4), and (iii) a single component with (0.7,0.7). Note that AiV(e) is 
normalized such that f AN(e)de = 1. 



4.2 Insights into the formation of the Milky Way 

In the previous subsection, we have demonstrated a non- 
constant radial /3-profile in the inner-halo component. In 
this subsection, we briefly discuss the insights of this result 
into the formation of the Milky Way halo. 

4-2.1 Radial /3 -profile as a probe of the relaxation process 

iLvnden-Belll l|l967f l argued that a stellar system that expe- 
rienced a rapid change of the gravitational potential would 
evolve into the relaxed system with an isotropic, ergodic dis- 
tribution function. He called this phenomenon 'violent relax- 
ation.' He expected that if the relaxation process is spatially 
limited in the central part of the system, the resultant distri- 
bution function would be altered so that it is isotropic only 
at the center, while staying anisotropic in the outer part of 
the system. As an example that could approximately repre- 
sent such a system, he mentioned a distribution function of 
the form f(E, L) — f(Q), which was later sh own to accom- 
modate a rad ial profile of fit — r 2 /(r 2 + r 2 ) l|Osipkovlll979l ; 
lMerrittj|l985T ). In his view, the scale radius r a corresponds 
to the 'relaxation radius' inside which the relaxation process 
acts effectively, in reasonable agreemen t with the numeri cal 
experiments of violent relaxation (e.g.. Ivan Alba da 1982). 

The distribution function we used in this paper [equa- 
tion (|12[) ] is a natural generalization of f(Q), and allows a 
free parameter of /3o at the center of the system in addi- 
tion to r a - If the actual /3-profile is described by the form 
in equation (|13[) . estimates of r a or /3o may be used as a 
guide to understand the relaxation process that has acted 
on the system. For example, if the scale radius r a could in- 
dicate the spatial r each of effective relaxation, as noted by 
ILvnden-Belll (|l967h . our present constraint of r a < 10 kpc 
or equivalently /3o < 0.5 for the inner halo (section I4.1.3[) 
might serve as a hint to uncover a trigger of such relaxation 
in the early Galaxy. Moreover, if some notable deviation of 
/3o from zero could be seen, it might indicate that the re- 
laxation would not be effective enough to completely erase 



the initial dynamical condition even at the center of the sys- 
tem. In order to impose a stronger constraint on j3o than our 
present constraint of /3o < 0.5 (sectionl4.1.3"|). it is necessary 
to re-analyse a sample by ICarollo et alj (|2010t ) and derive 
AA r b s (e) in a spherical halo potential, with which we can 
directly compare our models. 

We note that our present constraint on the large-scale 
profile of ft beyond the solar neighbourhood should be tested 
against direct measurements of j3 at various radial distances 
away from the center. In fact, I Bond et al.l (2010) reported 
that such analysis for a sample of stars with [Fe/H]< —1.1 
in 3 kpc < R < 13 kpc and 1 kpc < \z\ < 5 kpc supports a 
nearly constant profile of f3. More works for direct measure- 
ments of /3 are obviously necessary for further discussions 
on its radial profile in the inner halo. In the meanwhile, we 
plan to extend our present analysis to include another type 
of distribution functions in which /3q > /3q is also allowed 
(Hattori 2011, in preparation). With these calculations, we 
would be able to t est whether /3(r) is a decreas ing function of 
r (as proposed bv lSommer-Larsen et aljfl 997f ) or not. How- 
ever, we stress that our present calculations of AJV(e) do 
help constrain the large-scale profile of ft from a sample of 
solar-neighbour halo stars. 



5 CONCLUSION 

In this paper, we formulate the eccentricity distribution 
of the solar- neighbour halo stars AiV(e), with simple as- 
sumption of the spherical halo. By adopting two appropri- 
ate halo potentials, we show that AJV(e) is highly depen- 
dent on the radial profile of velocity anisotropy parameter 
/3. Moreover, we show that our theoretical calculations of 
Ai V(e) are usefu l in exp laining some properties of AiV b s (e) 
bv lCarollo et al.l |2010l ). such as the linear e-distribution in 
the low-e region and the metallicity-dependent fraction of 
high-e stars. 

We have demonstrated that the observ ed fraction of 
high-e stars with -2.0 <[Fe/H]< -1.5 by ICarollo et all 
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is smaller than that expected from our models with 
a constant profile of /3 in the inner halo. This result places a 
constraint on its radial profile for the inner-halo component 
of the Milky Way, yielding /3o < 0.5 at the galaxy center in 
contrast to the obsered value of /?© ~ 0.7 in the solar neigh- 
bourhood. This result further shows that the scale radius 
r a of the /3-profile should be smaller than ~ 10 kpc, which 
might imply that some relaxation process that acted on the 
inner-halo component was effective only within this radius 
away from the center of the Milky Way. 



ACKNOWLEDGEMENTS 

We thank Beers, T., Carollo, D., Minezaki, T., Tsujimoto, 
T., Yamagata, T., Sakata, Y., Kakehata, T., and Fujii, H. 
for useful discussions and suggestions. KH is supported by 
JSPS Research Fellowship for Young Scientists (23-954), and 
partly by Hayakawa Sachio Foundation. 



REFERENCES 

Abadi, M. G., Navarro, J. F., & Steinmetz, M. 2006, MN- 

RAS, 365, 747 
Bond, N. A., et al. 2010, ApJ, 716, 1 
Brown, W. R., Geller, M. J., Kenyon, S. J., & Diaferio, A. 

2010, AJ, 139, 59 
Carollo, D., et al. 2010, ApJ, 712, 692 
Chiba, M., & Beers, T. C. 2000, AJ, 119, 2843 
Chiba, M., & Yoshii, Y. 1997, ApJ, 490, L73 
Chiba, M., & Yoshii, Y. 1998, AJ, 115, 168 
Cuddeford, P. 1991, MNRAS, 253, 414 
Eggen, O. J., Lynden-Bell, D., & Sandage, A. R. 1962, ApJ, 

136, 748 

Gilmore, G, Wyse, R. F. G, & Kuijken, K. 1989, ARA&A, 
27, 555 

Hattori, K., & Yoshii, Y. 2010, MNRAS, 408, 2137 (HY) 
Henon, M. 1959, Annales dAstrophysique, 22, 126 
Lynden-Bell, D. 1960, MNRAS, 120, 204 
Lynden-Bell, D. 1962, MNRAS, 124, 1 
Lynden-Bell, D. 1967, MNRAS, 136, 101 
May, A., & Binney, J. 1986, MNRAS, 221, 857 
Merritt, D. 1985, AJ, 90, 1027 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 
490, 493 

Norris, J., Bessell, M. S., & Pickles, A. J. 1985, ApJS, 58, 
463 

Osipkov, L. P. 1979, Soviet Astronomy Letters, 5, 42 
Sommer-Larsen, J., Beers, T. C, Flynn, C, Wilhelm, R., 

& Christensen, P. R. 1997, ApJ, 481, 775 
Sommer-Larsen, J., & Zhen, C. 1990, MNRAS, 242, 10 
van Albada, T. S. 1982, MNRAS, 201, 939 
Voglis, N. 1994, MNRAS, 267, 379 
Yoshii, Y., & Saio, H. 1979, PASJ, 31, 339 
Zhao, D. H., Jing, Y. P., Mo, H. J., Borner, G. 2003, ApJ, 

597, L9 



APPENDIX A: LINDBLAD DIAGRAM 

The aim of this appendix is to show the phase-space region 
in which observable stars are distributed and to demon- 



strate how the cutoff eccentricity e cu t depends on specific 
energy E. 

Let us consider a spherical system of stars in a spherical 
gravitational potential V(r), and an observer 0\ whose 
survey region is r_ < r < r+ As noted in section 
12.31 observable stars for Oi are defined as those satisfying 



< r a 



and 



< r+. Since 



and 



depend 



on E and L, the phase space region which is occupied by 
observable stars can be clearly sh own in the (E, L)-ph ase 
space, called Lindblad diagram ' (|Mav fc Binnevlll986l ). 
Figure lATl is an example of such diagram for the Keplerian 
potential. We see from this diagram that stars with E < 
V(r~) are unobservable (AT r — 0) regardless of e (e cut = 
1) because r apo < r_, while stars with E c (r~) < E < 
E c (r+) are observable (AT r > 0) regardless of e (e cu t = 
0) because both r_ < r apo and r pcr i < r+ are satisfied. 
On the other hand, stars with V(r ) < E < E c (r~) are 
observable only when e > e cut where e cut corresponds to 
r apo = r_ . In this range of E, e cut decreases monotonically 
from 1 at E = V(r-) to at E — _E c (r_). Similarly, stars 
with E c (r+) < E < are observable only when e > e cut 
where e cut corresponds to r pcr i = r+. In this range of E, 
e cut increases monotonically from at E — E c (r+) to 1 
at E = 0. 



12 This assumption of survey region may seem too simple, but is 
justified by considering a more realistic observer O2 at r = (r+ + 
r_)/2 who observes stars within a distance (r+ — r_)/2 away 
from him. Since the survey region of observer Oi includes that of 
observer O2, the latter region well represents the former as long 
as (r_)_ — r_)/2 is small. Accordingly, the eccentricity distribution 
for these observers would be similar to each other in this case. 



Halo stars in the solar neighbourhood 




E / (GM/r ) 

Figure Al. Lindblad diagram with the Keplerian potential 
V(r) = —GM/r for an observer at r = tq whose survey re- 
gion is r_ < r < r+. In this case, we assume r_ = O.8r0 and 
r_|_ = 1.2rQ. The abscissa represents the specific energy E in 
units of GM/r@ and the ordinate represents the specific angular 
momentum L in units of ^/GMtq. No stars are allowed above 
the contour of e = (solid thin line), while observable stars are 
distributed below A-B-C-D (solid thick line), where A-B repre- 
sents the contour of r apo = T— and C-D represents the contour of 
r peri = r +- Here, the positions of A, B, C, and D correspond to 
E = y(r_), E c (r-), E c (r+), and 0, respectively. The contours 
of e = 0.3,0.6, and 0.9 are shown by dashed thin lines, which 
help understand that the cutoff eccentricity e cu t decreases with 
increasing E from A to B [V(r_) < E < E c (r— )], remains zero 
from B to C [£ c (r_) < E < E c (r+)], and increases from C to D 
[E c (r+) <E<0]. 



